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Process of two pion production in electron - polarized proton scattering is investi- 
gated. In Weizsacker- Williams approximation the differential spectral distributions 



and the spin-momentum correlations are considered. The spin correlation effects 
caused by p-meson widths are estimated to be of order of several percents. Both 
Qh' channels of it + it and it + it creation are considered. The effects of intermediate ex- 



cited baryons are not considered. The spectral distributions on pion energy fractions 
in polarized and unpolarized cases are presented analytically and numerically. 



> 

I. INTRODUCTION 

cn 

We consider the kinematical region of small angles electron scattering, which corresponds 
to the fragmentation region of the initial polarized proton. Besides we are interesting the case 
of small invariant masses of the final system consisting of the nucleon and the pair of mesons 
supposed to be not exceeding the mass of baryon resonances. In this case the intermediate 
state with p-meson becomes important. Apart of Weitzsaker- Williams (WW) enhancement 



-|3], which is associated with "large logarithmic factor" (In (s 2 / (M p 2 m^)) ~ 20), an ad- 
ditional enhancement associated with Breit-Wigner character of the intermediate p-meson 
propagator takes place when considering the total cross section as well as differential dis- 
tributions. The rather large width of p-meson induce the measurable one-spin asymmetries 
|| then initial proton is polarized transversally to the electron beam direction (laboratory 
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frame is implied). 

Our paper is organized as follows. We start the general consideration in Section HT1 where 
we consider the processes 

e~(pi)+p(p) -»■ e~(pi) +p(p') + vr + (g+) + 7r~(g_); 
e~(pi)+p(p) -»■ e - ^) + n(p') + vr + (g + ) + 7r°(g ). 

The definite channels of 7r + 7r~ production and 7r + 7r° production are considered in Sections lHIl 
[IV] correspondingly. In Section |V] we gather the final results and demonstrate some numerical 
estimations for experimental observables. 

Comparison with data obtained at HERMES facility are given. 



II. GENERAL CONSIDERATION 



The elastic scattering kinematics is implied 



2 12 
Pi = Pl 



2 12 
P = P 



M 2 ; 



2 2 2 

q± = q =m . 



(1) 



Besides we consider the case of small angles of electron scattering: 



-q = -(Pi-Pi) < s 



s = 2pp 1 > M 1 



(2) 



So for HERMES facility we have s ~ 30 — 40 GeV 2 . Matrix element in Born approximation 
has a form 



M e P ->en2. = ^ { _ J u (q)g^, 



q 



(3) 



where J u (q) is the hadron current 



Mq) = u (p')°Mp), 



(4) 



which obeys the current conservation condition J\(q)q x = 0, operator O v will be defined in 
Sections fTTTl IIVI The quantity G is a product of p-meson coupling constants G = g p7Tn gpNN- 
For theperipheral kinematics the Sudakov parametrization of 4-vectors of the problem is 
useful p-l8[. For this aim we introduce two almost light-like 4 vectors pi and p: 



Pi=Pi- P- 



p\ = O [m 



S 

,M 2 



M 2 



P = P~Pi- 



e o 
S z 



~2 



P 



O ( m 



M 4 



e o 
S Z 



i = s; 



2pp = M 2 



2piPi = rn e 



(5) 
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And the corresponding parametrization reads as 

q = a p + ppx + q±; 

q± = x±P + P±Pi + q±±; 

qo = %oP + 0oPi + qo±] 

p' = x p + fipx + q ± , (6) 

where the orthogonality conditions (pia± = pa± = 0) takes place. The 4- vectors a± have 
only spatial components, besides it's z-axis component is zero (we choice as a z-axis the 
direction of motion of initial electron a± = (0,0, a), where a is a 2-dimensional vector 
transversal to the z-axis. Thus we have 

a\ = -a 2 < 0. (7) 

The conservation law and final particles mass-shell conditions give the following relations: 

sfi± = — (£ + m 2 ) , 
x± 



x V 



J2 , A/r r 

X 

q = p* + q+ + q~, (8) 



X ~\~ X— ~h Xj^ — 1, 

s 

for the subprocess of tt + tt~ photoproduction: 7*(q r ) + p(p) —> p(p') + vr + + 7r _ . The similar 
relations for the subprocess 7* (5) + p{p) — > n(p') + 7r + + 7r° have a form 

s0 + = ^{f + + m 2 )- 
x + 

sPo = — (</o +1T1 2 )] 
Xq 

S p = I(^ + M 2 ); 



X 

q = P + q+ + %, (9) 

X + x + + x = 1, 
M 2 

+ — = 0+ + Po + 0'- 
s 

Significant simplification comes from the Gribov's form of the metric tensor g^ v in photon 
propagator jf|: 

gT = 9±+l Wp\ + p u pT\ ~ 2 -/pi- (10) 
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In the kinematics we are considering the only (2/ ' s)p^p\ component of this tensor is relevant. 
Other terms give the contributions which are suppressed by factors 

o(^)«, 

We also systematically omit the terms of order 

O (£) ~ °- 02 - (12) 



M 2 ) 

The terms of this order of magnitude (which are systematically omitted below) determine 
the accuracy of our consideration. 
Thus matrix element has a form: 

M e ^ eN2 * = *^N e Ml 2 = 7T + 7r-,7r + 7r°, (13) 
N e = -[u(pi)7 M «(Pi)]^- 



It is easy to see that 



J2\N e \ 2 = 2, (14) 



and that the quantity Nh does not depend on s in the large s limit. 
Let us now parameterize the phase volume of the final particles: 



jdVi d 3 p' d 3 qi d 3 q 2 
2E[ 2E> 2E 1 2E 2 



dT = (2n)-^^-^^5 4 (p + Pl -p'- pi - qi - q 2 ). (15) 



We use below the on-mass shell relation (i = +, — , 0) 

^J^ = d 4 q[5 (q 2 - m) = ^dxidotid 2 ^ (sa^gf - m 2 ) = -^—dxid 2 qi. (16) 
Introducing the additional unit factor 

1 = d 4 q S 4 (pi - p[ -q) = d/3 d 2 q 5 4 ( Pl - p[ - q) (17) 



we put the phase volume to the form 



dT = ( 27r )- 8 pl^d'qdY^qz. (18) 
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The channel of 7r + 7r~ production corresponds to q\ = <f_, Xi = x_, q 2 = q + , x 2 = x + ; while 
the channel of 7i + n° production corresponds to q\ = q , Xi = x , q 2 = q + , x 2 = x + . In phase 
volume transformation we used the mass shell conditions for scattered electron 

pf — ml = (pi — q) 2 — m 2 e = sa(3 — <f — sa — m 2 e (5 = 0, (19) 

and the conservation laws ([8]) and (J9j). Specific attention should be payed to expressing the 
transfer momentum squared q 2 . Using the scattered electron mass shell condition ffl9l) we 
can write it as 

g 2 = saP -? = -j^- (<f + m 2 J 2 ) « - (<f + ml(3 2 ) . (20) 

The quantity sj3 can be expressed in terms of invariant mass squared of final hadronic 
system: 

s = s$ = (p 1 + q, + q 2 f + <f- M 2 . (21) 

This quantity in Weitzsaker- Williams approximation can be written as 

M 2 

(So)q=0 = 1X1X2(1 -X) + V 2 X{\ - X) 2 + Xi(l - Xi)%2 + 

XX\X 2 

+ x 2 (l - x 2 )xi + 2xxx 2 x] = SlM , (22) 

XX\X 2 



where 



m q 2 q x q 2 



M' AJ M 2 ' M 2 ' 

We see that the transferred momentum q is the space-like 4-vector with non-zero module. 
One can see that the module of the matrix element tends to zero in the limit then transferred 
momentum q goes to zero. To be convinced let us write the quantity N h as 

M h = -PiUq) = ~(q - q±y-Uq) = --<£•/„(?), (23) 

S So So 

where we have used the current conservation condition. 

Here and further we restrict ourselves only by the approximation of a "large logarithm" 
(i.e. so called Weitzsaker- Williams approximation). For this aim we perform the integration 
over the transferred momentum using the relation 

d 2 q Mi 1 



7r(<f + m 2 J 2 ) 2 2 



<%(£-!) 
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where L is the " large logarithm" : 

2 s 2 s 2 

L = ^^l ^23. (24) 

Writing the matrix element of the pions photoproduction sub-process j*(q) +p — > Nttitt2 

as 

M l h = —u(p')diu(p), i = n + n-,n + 7i°, (25) 
So 

and using the density matrix of a polarized proton as u(p,a)u(p,a) = (p + M)(l — 750), 
where 4- vector a is the polarization vector of the initial proton (i.e. p ■ a = 0), we obtain for 
the differential cross section: 

^ M " d ^ _ <*& (L _ 1)xx , £ (26) 

dx 1 dx 2 d 2 q l d 2 q 2 6AMV [ L ) xx ^ ^ W 

where 

1 



Si = 4i^ Tr 



{p' + M) Oi (p + M) (1 — 75a) Oi . (27) 



We note that for HERMES conditions the factor a = a 2 G 2 (L - l)/(64M 2 7r 3 ) have a rather 
large value o"o ~ 2 ■ 10~ 30 cm 2 , i.e. of order of two microbarns. The explicit expressions for 
the quantity Si for both processes considered here are given in Appendices [A], [Bl 

III. MATRIX ELEMENT OF vr+vr PRODUCTION ON A PROTON 

The matrix element of tt + 7t_ mesons production can be presented as a sum of two sepa- 
rately gauge-invariant contributions 

M n+7T _ = M 1 + M 2l (28) 

First term describes the subprocess of creation of pion pair by (virtual) photon and (virtual) 
p-meson 

x r ((2g_-g),e)(-2g + + gi) p + (2g_ - qi ) p {{-2q + + q) , e) _ ^ ) 
I {q — q_) 2 — m 2 {q — q + ) 2 — m 2 p ) ' 
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where e = e(q) is the polarization vector of virtual photon, qi = q + + g_ — q, q 2 = q + + g_. 
Second term describes the subprocess of Compton scattering on a proton with emission of 
p-meson which subsequently decays into pion pair 
G 



Mo 



q! - Ml 



, p'-q + M p + q + M , , , 
U \P ) i e T~i ~ ^ v + v 7~~< ~ T7b-. e ( U (P) 



(30) 



(p'-q) 2 -M 2 (p + q) 2 -M 2 

where v = —q+ + q-- Here we omit the QED coupling a since it is already taken into account 
in flU. 

Both contributions (Miand M2) satisfy the gauge condition: they turn to zero when use 
the mass shell conditions and replace e(q) — > q. 

The matrix element (1251) can be written in the form 

M v+V _ = — \u{p')6 + _u{p)\ , (31) 

where we left only quadratic over q terms and operator + _ has a form: 



A 



Aq + + Bq_ + Cq± + Dv + Eq^piV + Fvpxqx, 



where 



A = 



4x 2 x + qq_ 
Sl D p M 2 ' 



B 



4x 2 x_ qq + 
Sl D p M 2 '' 



D = 



2x 2 _x 2 + fq_ 
Sl D p M 2 ' 



2x 



X + X- 



and the p-meson and pion propagators denominators are 

D p = xfi 2 + z/ 2 (l - xf + X1 + X2 + 2 X , 

= -/u 2 (l - i7 p )x + x_ + x\x2 + x\xi - 2x+x-x + v 2 {l - x) 2 , 



(32) 



7p 



M 



rn 
M' 



IV. MATRIX ELEMENT OF vr+vr PRODUCTION ON A PROTON 



For the subprocess 7*(g) +v(p) ~~ >* n (p') 7l +( ( l+)' n 'o{ ( lo) the amplitude M +0 contains besides 
the Compton-like amplitude 



Mi 



G 



(q + + q ) 2 -M 2 



u(p')(q - 9h 



p + g + M 
(p + g) 2 - M 2 



eu{p) 



(33) 
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also three terms corresponding to conversion of virtual photon and virtual p-meson to pion 
pair: 

u(p'hp u (p) 



Mo = G 



x 



M 3 = -G 
M A = G 



((<?+ + go) 2 - M p 2 ) ((g+ + q - g) 2 - M 2 ) 
ev+ (~q+ - go - q) p + v+p (2q + + 2q -q)e + e p v + (2q - q+ - g ) 
u(p')(-q+ + % + q)u(p)e(-2q + + q) _ 



((g + + g ) 2 - Mf) - g) 2 - m 2 ) 
u(p')eu(p) 



(g + + g -g) 2 -M2- 
It can be checked that the sum 



v = v + = -q + + g . 



M+o = Mi + M 2 + M 3 + M 4 , 



(34) 



obeys the gauge condition. Writing this amplitude in the form similar to (1311) leads us to 

sG 



where: 



O 



Mi 



1 



s 



u(p')O +0 u(p) 



(35) 



y+Piq± ^ 

sD n D P D^ 
1 



[-2fip x qv+ + 2v + q(q + + g ) + s 1 (x + - x )q±] 



[six + q± + 4gg + g ] 



(36) 



where the expressions for s\, D p , D n can be obtained from the relevant expressions for 
process of 7r + 7r_ production (see Section HTT1) by replacement g_,x_ — > qo,x . 



V. RESULTS 

The differential cross section of pion pair electroproduction processes which we are con- 
sidering here can be obtain using Equ. (|26|) . In case of 7r + 7r~ pair production we have to 
insert the quantity S* + _ into this formula (TSB"]) . This quantity has a form: 

s+ _ = (1 _ P±) j - 2 ^:S: ImD -A + iR t D - M l d \ + 
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x + 1 x_ 


0.2 


0.3 


0.4 


0.5 


0.6 


0.7 


0.2 


0.268 


0.359 


0.362 


0.312 


0.250 


0.138 


0.3 


0.337 


0.613 


0.708 


0.620 


0.410 




0.4 


0.351 


0.712 


0.899 


0.808 






0.5 


0.315 


0.625 


0.792 








0.6 


0.240 


0.406 










0.7 


0.137 













Tab. I: The function (x+,x-) (denned in (|4ip ) is presented for different values of final pions 
energy fractions x + and x_ for 7r + 7r~ production case. 

where permutation operator P± acts on any function F (x + , x_, q + , q_) in the following 
manner 



P±F (x+,x-,q+,q-) = F x+, <f_, q+) 



(38) 



and the quantities A, C, D, F in (1371) are given in Appendix lAl 

In case of 7r + 7r° pair production we have to insert the quantity S+o into this formula ( |26l ). 
This quantity has a form: 



S. 



2/mD, 



+o 



siD*\DJ* 



M V 1 2 p! M 



+ 



C 



2 :D + ^*E + 2 



D P \DJ 2 sMDJ* D P \DJ 



' D |D p G 



ID. 



(39) 



where quantities Al, A 2 , B, C, D, E, F, G, H are given in Appendix iBl 

Dalitz-distribution (the distribution over final pions energy fractions) can be written in 
the following form 

da a 2 G 2 



dx x dx 2 64M% 3 
where function _F unp (x 1; x-i) 



(L-l)F i unp (x u x 2 ) 



(40) 



^unp / x f d q\d q 2 Si 
F- l (x 1 ,x 2 )= / xxix 2 -~. 

J 7T 7T Si 



(41) 



This function in case of 7i + n production is presented in Table [J while in case of 
production is in Table [TTJ 



7T 7T 
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x+ 1 x 


0.2 


0.3 


0.4 


0.5 


0.6 


0.7 


0.2 


56.325 


23.089 


2.897 


0.418 


0.052 


0.003 


0.3 


23.166 


14.824 


5.237 


0.421 


0.016 




0.4 


2.876 


5.063 


0.990 


0.071 






0.5 


0.415 


0.419 


0.067 








0.6 


0.051 


0.016 










0.7 


0.003 













Tab. II: The function F"q P (x + ,xo) (denned in (|4ip) is presented for different values of final pions 
energy fractions x + and xo for ir + -K° production case. 

In the experimental setup when we fix the azimuthal angle if) between the proton po- 
larization vector a and transverse momentum of one of pions <f + , i.e. if) = (a, some 
polarization dependent contributions appear: 



d \q[\ 



d 2 q 2 



S; 



XX1X2-5 = \a\ simp Ff 01 (xi, x 2 ) + F" lip (xi, x 2 ) ■ 



71 



And thus in experiment the asymmetry 



A(x 1 ,x 2 ) 



F po1 (x 1 ,x 2 ) 



(42) 



(43) 



can be measured. This asymmetries for both channels are presented in Figs. [U |5] as a 
function x + for fixed values of x„ or xq. 

We should note that the meaning of quantities x, x±, x 2 is different in laboratory reference 
frame and in the center of mass of the initial particles (cm. frame). In cm. frame they 
are the energy fractions of the jet consisting of the recoil proton and pions which obey the 
conservation laws 



x + X\ + x 2 



2Ei 



2m 



< x, < 1, 



2M 



< x < 1. 



(44) 



In Laboratory frame, keeping in mind the explicit form of the light-like vectors p and 
Pi ~ Pi- 



M 

p=— (1,-1,0,0), 



Pl = E(l,l,0,Q), 



we have for energies: 




Fig. 2: The asymmetry A(x+,xo) (defined in (|43p ) as a function of x+ with xq fixed to definite 
values. 
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S i 771 

M+ = E 1 + E 2 + E', —<Xi<l. 

2Mxx 1 x 2 M 



The values of 3-momenta squares are 



2 -2 i 

= & + 



rA . 9 "l 2 



M of + m : 
-Xi + 



2 2Mxi 

2 i A#2\ 2 



P = (?i + 12) + — + 



2 2Mx 

We specifically emphasize that in left part of these equalities are 3-dimensional vectors and 
in the right part - the two dimensional vectors (their components perpendicular to the initial 
electron direction are zero). 

Dalitz-plot distributions for unpolarized case (see Fig. [T]and|2j while gq ~ 2 fib) seems to 
be accessible in experiments of HERMES. The one-spin asymmetries (see f )43|) ) have a value 
of several percents and as well can be measured. Preliminary comparison of these quantities 
are in qualitative agreement with data obtained at HERMES Jjj]. 

Appendix A: The production of 7r + 7r pair 

Here we present the the explicit expressions for the values A, C, D, F for process 7r + 7r_ 
production (see (1371)): 

A = A lX+ + -MX- + \Msi + -A 5 ; (Al) 

c = Cix+x- + \c-ix+x-x + c 4 x+ + \c & x 2 + \{C 7 + Sl C 8 )x+X- + 
(C 9 + Sl C w )x+X + (Ci 3 si + C 14 + C 15 sl) x+ + i(C 19 + C 2{j s\) x + -C 2l s\ + ^C 22 s\- (A2) 

L» = ^DiX+X- + D sX+X- + A>X+X + \d 7X +X-X + DsX 2 + + \x\Dio + SiD u ) + 
^X+X-Dn + (D 13 + sxD^x+X + {Dns x + D 18 + slD 19 )x+ + 

-(D 23Sl + D 24 + £> 25 s?)x + \d 2& s\ + -D 27 s?; (A3) 



F = + F 3X 2 + X- + F^X+X 2 + F 7 %+x + \f 9X +X~X + xl(F loSl + F n ) + 
± x 2 (F MSl + F l5 ) + ^ 



2 

1 X 2 (F 14 s 1 + F 15 ) + \x+X-{F 1&Sl + F 17 ) + (F 18Sl + F 19 )x+X + 



(F 22 si + F 23 + s\F 2A )x+ + \{F 2S + siF 29 )x + \s\F m + ^F 31 . (A4) 



The expressions for coefficients are 

Ai = —x(l — x)(l — x_)x z _; A 2 = xx + x_{l — x){2x + — x_(l — x + )); 
A 3 = 2x(2 — x-)x + x 2 _{l — x); A^ — x 2 _(l — x){l + x + — 

A 5 = xx 3 _x + (l — x)(x + — X-). 



d = -—x^x 2 , - x+(l - C 3 = 4x 4 (2 - x); C 4 = — xV; 

C 6 = 4x 3 ; C 7 = -2x 4 [x + x_(l - x) - 2[x\ + x 2 _)];C 8 = (1 - x) 2 ; 

x + x_ 

2 2 
C 9 = Ax 2 _x 4 (l + x-); Cio = x^x 3 ; C 13 = — x 2 _ x 3 (l - x); 

X-\. X-\- 

x 2 

C14 = 2x 4 _x A ; Ci5 = ; C 19 = Ax 2 _x 2 ,x 4 ; 

x + 

C 20 = —x 2 (l - x); C 21 = -x 2 (l - x); C 22 = ' 



Di = 2(xx^) x + ; D 3 = 2(xx-X + ) ; D 5 = x x_x + (l — x-) 
D 7 = x + x-x 2 (x- — x+)(2 — x); D% = x 2 x 4 _x + (l + xJ) 
D w = — 2{xx + xJ) 2 {l — x)(x- — x + ); D n = xx + x_{l + x)(x- — x + ) 

D 12 = -3(xx_x + ) 2 (x^ - x+) 
D 13 = x 2 _x 2 x + [x 3 _ + x 2 _(l - 2x + ) - x_x + (l + x + ) + 2x 2 + ] 
D u = -xx_[x 3 _ — x 2 _ + x_x + {<6 — x + ) + x + (3x + — 2)]; D 17 = —xx 3 _x + 

Z 

Dis = —x 2 x 2 + xt(x- —x + );Di 9 = ]-x_[x + — xxJ\ 

Z 

D 23 = -xx_x + (\ — x)(l — 3x)(x_ — x + ); D 2i = —(xx^x + ) 2 (x 2 _ + x 2 + )(x^ — x + ) 

Z 

D25 = ^x(l - x)(x- - x+); D 26 = ~x-x+(l - x)(x- - x+) 

D 27 = -x+). 
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F x = -xtx 2 + (l - F 3 = --x z _x\[x 2 _ — X-{\ — x + ) + x + (2 + x + )]; 

F 5 = 2xtx 3 + ; F 7 = x 2 + x 4 _(l + x + - F 9 = -2(x_x+) 3 ; 
F 10 = ~\xtx\{2 - F u = -^rix 3 [x_(3 - x) - 2z+]; F 14 = -2(x_a: + ) 3 ; 

Fib = 2(x_x + ) 4 (l - x); F 16 = i(x_x+) 2 [x 2 _ + x\ + 2(1 - ar)]; 

Fyj = -i(x_x + ) 3 [(l - x)(x 2 + x\) + 2(x_ - x + ) 2 ];F 18 = xtx+(x- - x + ); 

Fig = -xix 3 (3 - - x + ); F 22 = (x_x=) 3 (x_ - x+); F 23 = i(s_x + ) 4 (a;_ - x+); 

F 24 = — x^x+l— x 2 _ (1 — x_) — x+(3 — x_) + 2x + (l + x 2 )]; F 28 = (x_x + ) 4 (x_ — x+) 2 ; 

F 29 = --^(x_x + ) 2 [l + x 2 }; F 30 = --^(x_x+) 2 (l - x)(x- - x + ) 2 ; 
2x Ax 

F 31 = ^-x„x + {x_ - a; + ) 2 (A8) 



Appendix B: The production of 7r + 7r° pair 

Here we present the the explicit expressions for scalar coefficients in trace for the process 
7r + 7ro production (see (J39l) ): 

M = si\x+a>i + XoG2 + xa 3 ] + s\a^\ 
A 2 = x+ a 5 + Xoae + X a i + sia$ + a 9 , (Bl) 



with 



with 



1 1 

a\ = x 2 Xq; a 2 = x + x x 2 ; a 3 = 4(x x) 2 ; a 4 = — xx^; 

05 = x 2 Xq; a6 = XqX + x 2 (x + — Xq + x + xq); a 7 = —2x + xlx 2 , 

a s = x(x + — Xq); a 9 = x + x^x 2 (x + — x ). (B2) 

fi = s 1 (6 lX + + & 2 Xo + &3X) + &4S 2 , (B3) 

bi = —(xox) 2 ; b 2 = (xx + ) 2 ; 6 3 = —Ax + xqx 2 ] 64 = — xx + . (B4) 

C = C1X+X0 + c 2 XoX+ + c 3 x+XoX + c iX+ + (o> + SiC 6 )x+Xo + (c 7 + Sic 8 )x+X + 

(C 9 S 2 + S1C10 + cn)x+ci 2 s 2 xo + (sic 13 + Ci 4 )s 1 x + 

+ S 2 Ci 5 + S 3 Ci 6 , (B5) 
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!; '° c 2 = —2(x+ - l)(x+ + x - l) 2 , c 3 = 4x (x + + x - l) 4 , 



X-\- 

c 4 = _ 2x o(^+ + ^)- 1 ) ^ Cs = _ 2 ^ x+ _ 2 )xl( x+ + Xq - l) 4 , 

X_j_ 

(x+ + x - l) 3 (x+ + 2x ) 3f . 4 a; (x + + a; - l) 3 
6 = , c 7 = 4x^(a: + + x - 1) , c 8 = , 

X-j- X_|_ 

( *+ + X ° - 1)2 , c 10 = _ ^ + + ^o-l) 3 (^ + + 2xo) > cu = ^ + Xo _ 1)4? 



4a; + ' x + 

{X + + X -1) 2 (X++X -1) 3 

C12 = " , Cis = , C14 = X+X (X+ + X - 1) , 

ci 5 = + x - l) 2 (a;+ + x ), c i6 = - ^^^ — --(B6) 



£> = rfxx+ + c? 2 Xo + d 3 x 2 +Xo + d A xlx+ + 4x+X + 
4XoX + ^?X 2 X+ + 4x 2 Xo + d 9 x + XoX + 

(sidlO + rfll)x+ + (Sldl2 + dl3)Xo + (Ml4 + ^15)X 2 + 

(sidie + rfi 7 )x+Xo + (Mis + dig)x+X + (M20 + d 2 i)xoX + 

(S 2 (i 2 2 + M23 + ^24)X+ + (S^25 + M 26 + d 27 )X0 + 

(s 2 d 2 s + d w )x+sl d 3Q + sld 31 , (B7) 



16 



with 



1 1 

di = --x 2 xl(x - 1), d 2 = --x 2 x 2 + (x + - 1) 

d 3 = --x 2 x (xl + (x + - l)x + x + {x + + 2)) 
d 4 = -\-x 2 x + (x 2 Q + (x+ + 2)x + x + (x + - 1)) 

Z 

d 5 = x 2 xl(x + — x + 1)xq, d 6 = — x 2 x 2 + {x + — x — 1) 
d 7 = 2x 2 x + xl, d 8 = 2x 2 x 2 h xo, dg = —2x 2 x + xo, d w = -xx (x + + 2x — 2) 

Z 

1 1 
du = --x 2 x + xl(x + (x - 2) + x (x + 2)), d 12 = -xx + (2x + + x - 2) 

Z Z 

d 13 = -\ > x 2 x 2 + xq(x + (x + + 2) + x (x - 2)), du = -2xx + x , d lb = 2x 2 x\x\{\ - x) 
z 

d\% = — x((x+ — l)x — X + ) 

d 17 = -^-x 2 x + xq(xI + (x+ + 2)xl + (x+ - A)x + xq + x 2 + (x + + 2)) 
z 

d 18 = -xx (x + - x ), d 19 = x 2 xlx + (x + - x )(3 - x) 

d 20 = xx + (x + - x ), d 21 = -x 2 x 2 + x (x + - x )(3 - x) 

I f x 2 \ 1 

d 22 = ^ ~ 3x + + 2 J ' rf 23 = --xx + x (x + - x )(2 + x) 

1 1 / x 2 

d 2A = -x 2 x 2 + xl(x + - x ) 2 , d 25 = - [~~ - 3a; + 2 

1 1 

d 2 e = ~-x + xx (x + -x )(x-l), d 27 = -x 2 x 2 + xl(x + - x f 

d 28 = -x, d 29 = x 2 x 2 + (x + - x ) 2 , d 30 = - - x) 

= (j^i 

Axx + xq 



E = eix+Xo + e 2 xlx+ + e 3 xlx + e^Xo + e 5 X+XoX + e 6 X+ 
(sie 7 + e 8 )x 2 + (sie 9 + e 10 )xoX+ + (sien + e 12 )xx+ + 
(siei 3 + ei 4 )xXo + {s\e lb + siei 6 + ei 7 )x+ + (siei 8 + ei 9 )siXo + 

{s\e 20 + s x e 2 i + e 22 )x + s\e 23 + s?e 2 4, (B9) 
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with 



e\ = —xox 3 (xo + 1), e 2 = —x + x 3 (x + — 1), e 3 = — x + x 3 (x + — 1), 
e 4 = — 2x 3 x + x , e 5 = — xx (2x + — x — 1), e 6 = x 3 x 3 j (x + 1), 

e 7 = 2x 2 xo, es = — 2x 3 x + xl, eg = — x 2 , eio = xx + x$(x + + (x + — 2)xo), 

xox 2 (2x + — x ) n o 1 ^ . 

en = — , e 12 = -x x {2x + - x - 1), e 13 = --x {3x + + 2x - 2), 

e u = x + x x 3 (x + + (x+ - 2)x Q ), e 15 = — — , e 16 = -\x 2 Q x 2 (x Q + 2), 

Ax + 2 

3 3/ \ X{X + — Xq) 1 2 2 

e i7 = x + a; a; (x+ - x ), ei 8 = , e 19 = -x x,x , 

Axq 2 

e 20 = — ^2_ ; e21 = x 2 x (x+ - Xq + x ), e 22 = x 3 x+Xo(a; + - x ), 
2x-\- 

e 2 3 = \x(x + - x )(l - x), e 24 = j (— -J . (BIO) 



^ = fiX 2 + + /2X0 + fsX 2 + faxXo + fbX+Xo + faXoX 

(si/7 + /s)x+ + (Si/g + /io)xo + six/11 + S?/i 2 , (Bll) 

with 

/1 = -^x+xgxfco - 1), /2 = -^^x x(x+ - 1), / 3 = 2x\x\x, 

= — — X+XoX^X 2 ^ + Xq + 1 — x), /s = X+XXq^X-)- — Xq), fg = — X^XqX^X^. — Xq), 

h = 7^0(^0 + ( x + ~ 1 ) x +), fs = ^x 2 + x 2 x(x + - x ), f 9 = ^x + (x 2 + + (x - l)x ), 

f W = -^X%X%x(x+ - Xq), f U = X + X Q X, f 12 = ^XqX. (B12) 



with 



G = 9iX 2 + 92X0X+ + 9sXX+ + 94XXo + 
(si9 5 + 9e)x+ + SiXo07 + (sxgs + g 9 )x + s\g 1Q , (B13) 



111 1 

g x = -x + x 2 x 2 , g 2 = -x + x x 2 , g 3 = -x 3 Q x 2 , g A = -(x+ - l)x + x x 2 , g 5 = -x%x, 



1 : 



g 6 = -x + xlx 2 , g 7 = --x 2 + x, g 8 = -(x+ - x )x x, g 9 = --x + x 3 x 2 , g 10 = -(x+ - x jBU) 
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H = hx+ + h 2 Xo + hx + h 4 si + h 5 , 



(B15) 



with 



111 

--x+xl(xl + x + (x - 1)), h 2 = --x 2 + x (x 2 + + (x+ - l)x ), h 3 = --x 2 + xlx, 

X + (X + — X ) 2 X 1 2 2 2 . . 

= Ax ' = 4 X + X °( X+ ~ x °> ' ( * 
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